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ABSTRACT 

We investigate the nucleosynthesis in the neutrino-driven winds blown off from a SMq massive 
proto-neutron star (mPNS) temporarily formed during the collapse of a IOOMq star. Such mPNSs 
would be formed in hypernovae. We construct steady and spherically symmetric wind solutions. 
We set large neutrino luminosities of ^ 10®^ erg s“^ and average energies of electron neutrinos and 
antineutrinos in the ranges of £[/,,= 9 — 16 MeV and = 11 — 18 MeV based on a recent numerical 
relativity simulation. The wind solutions indicate much shorter temperature-decrease timescale than 
that of the winds from ordinary PNSs and, depending on tn, the winds can be both neutron-rich 
and proton-rich. In the neutron-rich wind, the r-process occurs and the abundance distribution of 
a fiducial wind model of the mPNS gives an approximate agreement with the abundance pattern of 
metal-poor weak r star HD 122563, although the third-peak elements are produced only when the Ve 
energy is much larger than the Ve energy. In the proton-rich wind, the strong i/p-process occurs and 
A > 100 nuclides are synthesized. The synthesized nuclei can be neutron-rich in some cases because 
the large neutrino luminosity of the mPNS supplies a sufficient amount of neutrons. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - supernovae: general - stars: abun¬ 
dances - stars: neutron 


1. INTRODUCTION 

The neutrino-driven wind from nascent neutron stars 
in core-collapse supernovae (SNe) is a possible site of 
heavy-element synthesis. A number of studies about 
the r-process, which produces about a half of elements 
heavier than iron, in the neutrino-driven wind have been 


made in the past decades (e.g., 

Oian WooslevI 119961: 

Hoffman et al.lll997tlOtsuki et al. 

200(1: iThomnson et al.l 

20011 IWanaio et al.l 120011: IRoberts et al.l 201011. How- 

ever, it was revealed that the 

strong r-process that 


makes all nuclei up to the third-peak elements {A 
195) hardly occurs in the standard condition of the 
neutrino-driven wind. The most critical issue toward 
the successful r-process is that the electron fraction Ye 
of the neutrino-driven wind is larger than 0.5, i.e., the 
wind is not neutron-rich as recent numerical simulations 
of proto-neutron star (PNS) evolution suggested (e.g., 
iPischer et al.ll2010l : lHudeDohl et al.ll2010ll . Electron frac¬ 
tion in the neutrino-driven wind still has uncertainties 
due to the input phy sics such as mean field effects on neu¬ 
trino opacities le.g.. iHorowitz et 'mi2012HRoberts et al.l 
[Mil . so the r-process might occur in the wind. Never¬ 
theless, recent studies have indicated that the r-process 
in such a win d is too weak to produce the third-peak 
elements le.g.. IWanaioll201,ll . 

Instead, the proton-rich wind is a promising site 
of the i^p-process, which is the nucleosynthesis pro¬ 
cess found in the 2000’s. Recent studies show that 
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some light proton-rich nuclei can be made in hot bub- 
bles and neutrino-dr i ven winds through the z^p-process 
(lErohlich et al.l 120061 : iPruet et'al] 120061 : iWanaiol 120061 : 
iWanaio et al.ll20lHl . 

In the 1990’s, hypernovae (HNe), more energetic 
explosive events than ordinary SNe, w ere discovered 
(iGalama et al.lll9^ llwamoto et ahlfTO^ . Although the 
explosion mechanism of HNe is yet unclear, there are 
some observational hints about their progenitors. First, 
progenitors of HNe are thought to be massive stars, with 
the initial mass larg er than ^ 25Mq (e.g., HN branch in 
iNomoto et al.l[^06l l. which is larger than the masses of 
the progenitors of ordinary SNe. Second, at least some of 
HNe s how the features of bipolar explosion (jMaeda et al.l 
l2006alibll . which indicates that the rotation of the pro¬ 
genitor may play a role in the explosion. Based on these 
studies, we assume here that HNe are bipolar explosions 
associated with the death of very massive stars. 

After the explosion, the neutrino-driven wind may be 
launched toward the direction of the bipolar explosion 
if a PNS is formed at least for a while. Indeed, ac¬ 
cording to recent numerical simu lations of the collapse 
of rotating massive stell ar cores ([Sekiguchi et al.l 120121 : 
iCerda-Duran et iini2013l l. a massive proto-neutron star 
(mP NS) is temporarily formed after the collapse (see 
also iSumivoshi et 1111120071120081 for spherical collapse). 
The mPNS is supported against the collapse by its ther¬ 
mal pressure and rotation and eventually collapses into 
a black hole. The lifetime of the mPNS is a few sec¬ 
onds even in the case of the collapse of a IOOMq star 
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l|Sekiguchi et alJ[2012D for an equation of state (EOS) 
with which the maximum mass of the neutron star is 
larger than 2Mq. For a less massive progenitor, the life¬ 
time would be much longer. 

Based on the above picture for HNe and the collapse 
of a massive star, we consider neutrino-driven wind from 
the mPNS in the HNe. We propose the HN wind model 
as follows: We assume that an explosion occurs toward 
the polar region when a massive star collapses. After the 
explosion, the neutrino-driven wind is launched from the 
mPNS surface toward the po lar region. For the mP NS 
model, we adopt the result in iSekiguchi et al.l (j2012f l. 

There are two main differences between HNe and or¬ 
dinary SNe. One of them is the (temporal) formation 
of more massive PNSs in HNe for the consequence of 
the collapse of the massive star. The other is the higher 
neutrino luminosity in HNe, which is expected to be pro¬ 
vided by the stellar-gas accretion to the mPNS along the 
directions that are not the direction of the bipolar explo¬ 
sion (e.g., equator in rotational coll apse). In addition to 
these, as some models of HNe le.g.. iBurrows et ffll2007tl 
suggest, the rotation and the magnetic field of the PNS 
could be important for HN explosion. We will discuss 
these effects in Section [5l 

In this work, based on this scenario, we investigate 
the nucleosynthesis in the HN wind. It is known that 
the entropy of the neutrino-driven wind strongly de¬ 
pends on the mas s of a PNS because of the general 
relati vistic effects (jCardall fc Fuller! Il997t lOtsuki et al.l 
12000(1 . Generally speaking, a high-entropy wind re¬ 
sults in more efficient heavy-element synthesis because 
less seed elements are produced through triple-a an d 
^He(an, 7 )®Be(a,n)^^C reactions (iHoffman et aLlHO^ . 
When the abundance of the seed elements is small, 
the neutron-to-seed ratio becomes large and thus fur¬ 
ther heavy-element synthesis proceeds. The expansion 
timescale has a crucial role on heavy-element synthesis. 
A shorter expansion timescale also causes less effective 
seed element production. According to the numerical 
simulation of the collapse of the very ma ssive star and the 
mPNS formation (|Sekiguchi et al.lF20I2ll . the neutrino lu¬ 
minosity of the mPNS is much higher than that of an 
ordinary PNS as mentioned above. In this case, the ex¬ 
pansion velocity of the neutrino-driven wind of the mPNS 
is also expected to be faster than that of the PNS, so the 
timescale of seed element production is also shorter than 
that of the PNS. Thus, the neutrino-driven wind from an 
mPNS, although the lifetime of the mPNS may be short, 
is expected to be an interesting site of heavy-element 
synthesis. 

In this work, we investigate possible nucleosynthesis 
processes to produce heavy elements in the neutrino- 
driven winds from mPNSs. Our model of the neutrino- 
driven winds and the wind parameters are described in 
Section [21 A brief description of the nucleosynthesis cal¬ 
culation is given in Section [31 The results of the nucle¬ 
osynthesis calculation and details of the nucleosynthesis 
processes are given in Section H) Finally, discussions and 
conclusions are given in Sections [5] and El Throughout 
this paper, we use the unit in which the gravitational 
constant G, the speed of light c, and Planck constant ft 
are taken to be unity. 

2. NEUTRINO-DRIVEN WIND MODEL 


2.1. Basic Equations 

The basic equations of the fluid are the equations of 
momentum conservation, baryon number conservation, 
and energy conservation. We ignore the self-gravity of 
the fluid and the change of the mass of the neutron star 
due to the mass ejection through the wind and the gas 
accretion through the equatorial plane. For simplicity, 
we consider a steady and spherically symmetric flow, so 
that the gravitational fields can be treated as those in the 
Schwarzschild spacetime. We do not consider the effects 
of the rotation and the magnetic field of the PNS, which 
will be discussed in Section El Then, the equations can 
be written in the Schwarzschild coordinate as 


M = 4TTr^ pu, 

du du l + u^ — 2M/rdP M 
dr dr p{l + e) + P dr ’ 

ds f de P dp 

dr ^ \ dr p^ dr 



( 1 ) 

( 2 ) 

(3) 


where r is the radial coordinate, M is the mass of the 
PNS, M is the mass ejection rate, which is constant in 
the steady wind, u = dr/dr is the radial component of 
the four-velocity of the fluid element, r is the proper 
time of the fluid element, p is the rest-mass density, P 
is the pressure, T is the temperature, e is the specific 
internal energy, s is the entropy per unit mass, and q is 
the net heating rate per unit time per unit mass. These 
equ ations are the same as those in lOtsuki et al.l l)2000f) 
and iWanaio et all l)200Ih . In addition to these equations, 
we solve the evolution equation for the electron fraction 
Ye, which is fixed in their work. Ye changes owing to the 
interactions between nucleons and neutrinos: 

dYf. dYe 
dr dr 

— (^^i?p “t“ Ai/yi -t- Aec + Apc)Ye -l- (A L^n + Apc),(4) 

where App,Ay„,Aec, and Ape are the rates for the elec¬ 
tron (anti) neutrino absorption reaction by free protons 
(neutrons) and electron (positron) capture reaction by 
free protons (neutrons), respectively. These reactions are 
written below: 


lye+P^n + e'^, (5) 

jze-I-n ^ p-I-e“. (6) 

The details of the calculation of these rates will be given 
in Section 1221 

The net heating rate q appeared in Equation m 
is composed of the terms of heating and cooling in- 
volving five processes and can be decomposed as (e.g., 
iQian fc Wooslevlll996D 

Q — QeN T Qei' 9e+e“ Y QuVj (7) 

where qvN^ Qev, and q^i, are the heating rates due to 
neutrino absorption on free nucleons, neutrino-electron 
scattering, and neutrino-antineutrino pair annihilation 
into electron-positron pairs, and qeN and qe+e- are the 
cooling rates due to electron and positron capture by 
free nucleons and electron-positron pair annihilation into 
neutrino-antineutrino pairs. The details will be given in 
Section 12.21 































nucleosynthesis in neutrino-driven wind in HNe 


3 


We rewrite Equations O-® to obtain 
, 2 A 2 \ du de f 2A dP M 


-- 

r \ op p 


de , 2 A 2\dT -I 

(“ —=(7U- 


dT 


dr 


AdP 

~^df ■ 

Au dYe (dP de 
p dr \dT dYe 
AdP 


2u f de P 


\ Q'X ' ’ 


u dp 


^ E.'] M. 

' ^ dp p J ^ 


dYe 

dr 


de 

We 


-A 


r r 
dP 


dp 


A dP / ^ 

p dYe \ dp 


dYe 

— (App “t" \ifn “t” Aec “1“ Apc)Eg 


“t“ {^lAn A Ape), 

where A = A{r) is defined as 


(9) 


( 10 ) 


Air) = 


1 + — 2M/r 

l + e + P/p 


( 11 ) 


Here Cg is the sound velocity given by 


= 


- (^\ - ^ /^^VY- — 


Kp), ap*pdT\dT 




where the thermodynamic quantities e and P are thought 
to be the functions of p, T and Ye. 

To close the system of equations, we adopt an EOS by 
iTimmes fc: Swest-^ (|2000ll . which includes the pressure of 
nucleons, electrons, positrons, and photons. 

We calculate the proper time of the fluid element from 
the definition of the four-velocity as 


r dr' 

W = / -r-yr- 
J u(r') 


(13) 


We construct the wind solution by solving Equations 
©-Col from the gain radius i?gain, where the net heat¬ 
ing rate becomes zero (q = 0). Eor simplicity, we assume 
that the radius of the neutrinosphere R^, is equal to the 
gain radius i?gain- This assumption is justified by the fact 
that the difference between and i?gain is small com¬ 
pared with the scale of the system. We further assume 
Fe = 0 at the neutrinosphere. Thus, for fixed R^,po 
and the luminosities and average energies of individual 
neutrinos L ^,., e ^,., we obtain the boundary conditions at 
r = i?j/(= i?gain) for temperature and electron fraction at 
the neutrinosphere Tq, Yg.o such that the conditions q = 0 
and Ye = 0 are fulfilled simultaneously. The boundary 
condition for velocity uq is set to be such that the wind 
solution becomes transonic. 


2.2. Heating and Cooling Rates and Reaction Rates due 
to Weak Interactions 


The heating rate due to neutrino and antineutrino ab¬ 
sorptions by free nucleons is given by 

Qi/N ~ 9.84 Nx [(1 — ]y[gy -|- j^gy] 


p2 


MeV s"^ g“ 


(14) 


where TVa is the Avogadro constant, is the lumi¬ 

nosity of neutrino or antineutrino of flavor i in units 
of 10^^ erg s“^, _MeV is the average neutrino or an¬ 
tineutrino energy of flavor i in MeV defined by e^. = 
aJ{E/,.)/{E^.) {{E//.) denotes the nth moment of the en¬ 
ergy distribution of the neutrino of flavor z), and R^fi is 
the radius of the neutrinosphere in units of 10® cm. In 
deriving the above expression, we approximate that all 
radii of the neutrinospheres of individual neutrino fla¬ 
vors are the same. The factor [1 — gi{r)] represents the 
effect that we see the neutrinosphere in a finite solid an- 
fe considering the effect of bending neutrino trajectory 
Otsuki et al.ll200(1ll . and pi(r) is defined by 


9i{r) = \/l - 


rW l-2M/r 


r / I- 2M/R^' 


and $(r) is the redshift factor given by 


(j)(r) = 


'1 - 2M/R^ 

1 - 2M/r ■ 


(15) 


(16) 


The heating rate due to neutrino-electron scattering is 
given by 




;2.17 Na 

Ps 


^ 51 ]V[eY “1“ ^i-'e ,51 ,MeV “1“ y,51 ,MeV 


Ji^$(r)® MeV s-^g-i 


(17) 


*^1^,6 


where is the sum of the luminosity of all neutrino fla¬ 
vors except for electron neutrinos and electron antineu¬ 
trinos, TmbV is the temperature of the fluid in MeV, ps 
is the mass density in units of 10® g cm“®, and Cjy. is de¬ 
fined by e^^ = {E/. )/ {E^^) (here means the neutrinos 
except for electron neutrinos and electron antineutrinos). 
In the following, we choose the relation between and 
e^i as e^. Ri l-14ej. as in lQian fc Woosle^ (|1996ll . 

The heating rate due to neutrino-antineutrino pair an¬ 
nihilation into electron-positron pairs is given by 

Na 


(Cj/^jMeV T £j/e,Mev) T y ^MeV 


X ^^4>(r)9 MeV s-^g-^, 

PsRZfi 

where 52 is defined using gi as 

52(0 = [1 - giir)f[giW +4gi(r) -k5]. 


(18) 


(19) 


The cooling rate due to electron and positron cap¬ 
ture by free nucleons qeN and the cooling rate due 
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to electron-positron pair annihilation into neutrino- 
antine utrino pairs Qp+p- a re calculated using the expres¬ 
sion in [Fuller et al.fp985f) (hereafter FFN85). The cool¬ 
ing rate per unit time per unit mass due to electron and 
positron capture by free nucleons is given as 


In 2 fkBTY 

qeN=meNp,-—- - 

[ft] \me J 

(l-Ye)J 


-h 


1 - exp(7?^ - Cn - vf) 
YeJ 


pc 


1 - exp(77^ - Cn - vD 
1.96 X 10-2 MeV s"^ g'^ 
{l-Ye)J 


-h 


_ 1 - exp{r]^ - Cn- V^) 

Y,J 


pc 


1 - exp(r7j - Cn - ) 


( 20 ) 


where subscripts “ec” and “pc” denote electron capture 
and positron capture, respectively. Here rUe is the elec¬ 
tron mass, {ft) K, 1.02 x 10^ s is the /t-value of the tran¬ 
sition between neutrons and protons, r]^ = p-e/feeT and 
rj^ = are the chemical potentials of electrons 

(positrons) and electron neutrinos (electron antineutri¬ 
nos) divided by ksT in the case of electron (positron) 
capture by free protons (neutrons), and Cn is —lS./k^T in 
the case of electron capture by protons and is A/kBT in 
the case of positron capture by neutrons (A = nrin—mp « 
1.293 MeV is the mass difference between neutron and 
proton). J is so called the phase-space factor of this 
reaction rate and is defined in appendix [A] 

The cooling rate due to the pair annihilation of elec¬ 
trons and positrons into a neutrino -antineutrino pair 
is calc ulated using the expression in iCooperstein et all 
(jl986[i . According to their expression, the cooling rate is 
given as 


^e+e“ 


1 Gp{De -|- + Dt) 

p 97r5 


{kBT)^K 


Ri5.43 X lO-"^ 


MeV s"^ 


Ps 


^-1 


( 21 ) 


where Gp is the Fermi coupling constant, Ug, and 
Dr are given using the Weinberg angle 0w as Hg = 1 + 
4 sin^ 0w + 8 sin^ 9w, Dp = Dr = 1 —4 sin^ 0w + 8 sin'^ 9w 
and we use the value sin^ 9w = 0.23 in evaluating these 
expressions, and K is the phase-space factor of this cool¬ 
ing rate, which is also given in appendix lAl 
The reaction rates of electron and positron capture by 
free nucleons Age and Ape are given using the expression 
in FFN85. In their expression, the rates are given by 


Age — 


In 2 / ^bT 


{ft) \me J 1 - exp(r7^ - Cn - pi) 

I 


a.96x 10-2 


1 - exp{T]P - Cn- vH 


s-\ ( 22 ) 




In 2 fkBT 


1 - exp(77^ - Cn - VD 


Ril.96 X 10-2 


rp3 

^MeV 


I 

1 - exp(77,f - Cn 




s-\ (23) 


where {ft), rjff, pf, and Cn are dehned in (l20l) . and I 
is the phase-space factor of this reaction rate, which is 
defined in appendix!^ In the following, we set rjff =0. 

The rates of neutrino and antineutrino absorption by 
fr ee nucleons Xyn and Xpp are given using the expression 
in iQtsuki et al.l (|200ClD as 


App R:! 9.84 f Cp^^MeV — 2AMeV + 1-2—MAG^ 

\ epe.MeV/ 

,, Lp^,5i<^>{r)^[l - gi{r)] „_1 

^ n 2 ® > 

Xpn ~ 9.84 (ep^^MeV + 2AMeV + 1-2—MAL^ 

V epe.MeV/ 

„ Lp^, 5 i^{r)^[l - gi{r)] 

X -m s . 

2.3. The Parameters of the Wind 

The parameters of the winds are the luminosities 
and average energies of the neutrinos of each flavor 

the mass of the neutron star 
M, the radius of the neutrinosphere Rp, which is as¬ 
sumed to be the gain radius iigain, and the mass density 
at the gain radius po- 

Here we explain the fiducial parameter set and param¬ 
eter range adopted in this study. We set the parame¬ 
ters on the basis of the numerica l relativity simul a tion o f 
the collapse of a IOOMq star in iSekiguchi et al.l (1201211 . 
They showed that in the center of the collapsing star an 
mPNS with the gravitational mass of M Ri 3 Mq and the 
neutrinosphere with the radius of Rp « 15 km is tem¬ 
porarily formed. Within its lifetime, on order of seconds, 
the neutrino luminosity and the mass of the neutron star 
become nearly constant in time. In this phase, the neu¬ 
trino luminosity of the neutron star is Lp ^ 10 ®^ erg s-^, 
the ratio of the luminosity of the electron antineutrino 
to that of the electron neutrino is Lp^/Lp^ Ri 1 . 1 , and 
the ratio of the luminosities of p, and r neutrinos to the 
electron neutrino luminosity is Lp^jLp^ Ri 0.1. Thus, we 
consider the steady flow in this phase under these condi¬ 
tions M = 3 M 0 , Rp = 15 km, Lp^ = 1.5 x 10®^ erg s-^, 
Lp^/Lp^ = 1 . 1 , Lp^/Lp^ = 0.1 from the numerical rela¬ 
tivity simulation. 

The average energies of the electron neutrinos 
and electron antineutrinos tp^ are very important be¬ 
cause they determine the electron fraction Yg and, con¬ 
sequently, the yield of nucleosynthesis. These values 
have some uncertainty caused by the approximate treat¬ 
ment of the neutrino transport and the hnite den¬ 
sity e ffects on neutrino o pacities (e.g., iHorowitz et al.l 
120121 : [R oberts et al.l 2012f). In the numerical simula¬ 
tion (ISekiguchi et al.M 2012 ll . their average energies are 
Cp^ « 11 MeV and Cp^ ~ 16 MeV. Thus, we construct the 
wind solutions and calculate nucleosynthesis on trajecto¬ 
ries of the wind for the energy region = 9 — 16 MeV 
and cp^ = 11 — 18 MeV, respectively. We also set the 
hd ucial parameter s of and tp^ as 11 and 16 MeV as 
in ISekiguchi et all (|2012l) . The average energies of the 
other flavors of neutrinos, p and t neutrino and its an¬ 
tiparticles, Cp^, is less important in constructing the wind 


(24) 


(25) 
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solutions or nucleosynthesis results because the luminos¬ 
ity of them is much smaller than that of electron and 
electron antineutrinos. Thus, they are set to be, from 
the simulation, 25 MeV for all wind solutions. 

We set the mass density at the neutrinosphere po to 
be 10^^ g cm“^. This value is 10 times as large as the 
mass density in the neutrino-driven wind models of an 
ordina ry SN in lOtsuki et al.l (|2000ll and iWanaio et al.l 
(j2001tl . where it is set as po = 10^° g cm“^. This is due to 
the very large neutrino luminosity in our model, which is 
about 100 times larger than their models. At the surface 
of the PNS, dominant heating and cooling processes are 
neutrino absorption reactions by free nucleons and their 
inverse reactions dSD and dH), and the balance between 
heating and cooling gives the relation among the neutrino 
luminosity the average energy of the neutrinos 
and the average energy of the electrons and positrons 
captured by nucleons E^. as 

L^el oc El. (26) 

The average electron energy at the surface of the PNS Ef. 
is almost equal to the chemical potential of electrons pe 
because of the electron degeneracy. The average energy 
of neutrinos can be estimated as the surface temperature 
T. Using the relation pe oc for relativistic degener¬ 
ate electrons. Equation (EH) can be rewritten as 

oc oc p^T-^. (27) 

Assuming that the surface temperatures are almost the 
same between the ordinary PNS and the mPNS, the sur¬ 
face density of the mPNS is about 10 times larger than 
that of the ordinary PNS. 

Summarizing the above, we give the parameter range 
of the wind in this work, from the simulation, in Table 

[II 


2.4. Comparison of Wind Solutions 

Here we discuss the differences between our mPNS 
wind and the ordinary PNS wind. The mass of the mPNS 
in our model is ~ SMq, and its radius (the radius of the 
neutrinosphere) is about ss 15 km. On the other hand, 
the mass of the ordinary PNS is ~ 1.4M0 and its ra¬ 
dius is about w 10 km. Thus, the mPNS in our model 
has a compactness parameter {C = M/R k, 0.3) larger 
than that of the ordinary PNS {C ~ 0.2). It indicates 
that general relativistic effects become more important 
in the mPNS wind than that in the ordinary PNS wind. 
In addition, the neutrino luminosity in the mPNS model 
is about ^ 10®^ erg s“^, which is about 100 times larger 
than that in the ordinary PNS (about 10®^ erg s“^). 

We solve the wind equations, i.e.. Equations (I51)- (ITU1) . 
with the two parameter sets given in Table |TJ One is our 
fiducial parameter set of the mPNS model, and the other 
is the parameter set of the ordinary PNS (the average 
energ ies of neutrinos are determined based on iWanaiol 
I 2 OI 3 II . We discuss the differences in the entropy, the 
velocity, and the temperature between them. 

Eigure |T] shows the radial profiles of the entropy, the 
velocity, and the temperature for the wind trajectories of 
the mPNS wind and the ordinary PNS wind. The termi¬ 
nal entropy of the mPNS wind, which is about 130-140 
(here the term “entropy” means the entropy per baryon 
in units of fee), is almost the same as that of the ordinary 


PNS wind (see top panel of Eigure[T]). Qualitatively, this 
can be explained as follows. The change of the entropy 
per baryon, As, of the wind may be given as 


As oc 


dq AQ 


(28) 


where AQ = J dq is the total heat obtained by the fluid, 
and T~^ = (AQ)~^ f dq T~^ is the mean temperature 
weighted by q. AQ will be of the order of the grav¬ 
itational potential at the PNS surface ~ M/R^, = C. 
Therefore, the wind blown off from the mPNS, which 
has larger compactness parameter C, is expected to have 
larger entropy than that of the wind blown off from the 
ordinary PNS. However, the change of the entropy of the 
mPNS wind is suppressed by higher T than that of the 
ordinary PNS wind (see bottom panel of Eigure |T]) be¬ 
cause of higher neutrino luminosity. These two opposite 
effects cancel each other, and as a result, the terminal en¬ 
tropies of the mPNS and the ordinary PNS winds have 
similar values. 

The biggest difference between the properties of 
the mPNS wind and the ordinary PNS wind is the 
temperature-decrease timescale at T = 0.5 MeV, de¬ 
noted by TTfi. 5 - It is defined as 


TT,0.5 


L 


T=0.5/e MeV 


dr. 


T=0.5 MeV 


(29) 


where e is the b ase of natural loga rithm (this defini¬ 
tion is based on lOtsuki et al.l 12000(1 . When the tem¬ 
perature decreases below T Ri 0.5 MeV, free protons 
and neutrons are reassembled into alpha particles so 
tt. 0.5 represents the timescale in which the seed ele¬ 
ments (the elements heavier than ^^C) are produced from 
alpha particles. Therefore, the temperature-decrease 
timescale at the point is the useful quantity that repre¬ 
sents the inefficiency of the production of seed elements 
thro ugh triple-g and ^He (an. 7 )®Be(a, n)^^C reactions 
le.g.. [Boffman et "nillOOTtl . 

The bottom panel of Eigure [1] shows the temperature 
of the mPNS and the ordinary SN as a function of the ra¬ 
dius. We see that, in the ordinary PNS wind, the radius 
of the critical point, indicated by a diamond, is larger 
than the radius where the temperature drops to 0.5 MeV, 
i.e., the alpha-particle formation radius. On the other 
hand, in the mPNS wind, the relation of these radii is 
opposite and the temperature decreases to 0.5 MeV af¬ 
ter the velocity reaches its asymptotic value (see also 
the middle panel). Therefore, the temperature-decrease 
timescale at T = 0.5 MeV of the mPNS wind is ex¬ 
pected to be shorter than that of the ordinary PNS wind. 
The value of rT,o .5 of the mPNS wind is about 6 ms 
and indeed much shorter than that of the ordinary PNS 
(si 40—50 ms). This trend can be understood more quan¬ 
titatively as follows. We define the temperature-decrease 
timescale at any temperature tt as 


-1 


1 dT 
T dr 
1 dp 
3/9 dr 
I /2m 
3 


1 du\ 
u dr J 


(30) 
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Table 1 

The Range of Parameters, the Fiducial Parameter Set of the mPNS Wind, and Those of the Ordinary PNS Wind 


Model 

M 

(Me) 

R. 

(km) 

PO 

(g cm-®) 

b.e 

h.. 

(10®^erg s ^) 



(MeV) 


Parameter range 

3 

15 

10” 

150 

165 

15 

9-16 

11-18 

25 

Fiducial 

3 

15 

10” 

150 

165 

15 

11 

16 

25 

Ordinary PNS 

1.4 

10 

lOlO 

1 

1 

1 

12 

14 

14 


In this equation, we assume dT/dr < 0 and, from the 
first to the second line, we use the fact that the radiation 
component dominates most of the total entropy and that 
the entropy is almost constant in time in all regions where 
T < 10^° K (see top panel of Figure[T|). From the second 
to the third equation, we use the proper-time derivative 
of Equation (HI)- The value tt^.b represents the time 
during which the temperature decreases from 0.5 MeV 
to 0.5/e MeV, and thus it can be written as tt{T = 
0.5 MeV) approximately. As seen in the middle panel of 
Figure [U the velocity of the fluid at T = 0.5 MeV (open 
circle shown in the figure) is about u ~ 0.22c and the 
radius is about r k, 140 km in the mPNS wind, while 
these values of the ordinary PNS wind are u sa 2.4 x 
10“^c and r « 35 km. Therefore, the first term in the 
parentheses of Equation (l30|) of the mPNS wind, which 
is dominant in the parentheses, is about one order larger 
than that of the ordinary PNS. 

Fewer seed elements are produced as tt,o .5 becomes 
shorter, and then the production of larger mass num¬ 
ber elements in the subsequent process after the seed 
element production becomes more efficient. Thus, while 
the entropy is almost the same, the circumstance of the 
mPNS wind of short temperature-decrease timescale is 
more attractive for heavy-element synthesis than that of 
the ordinary PNS wind. 


3. NUCLEOSYNTHESIS CALCULATION 

The nucleosynthesis calculation for the neutrino-driven 
winds in the mPNS is started when the temperature of 
the fluid element decreases to Tg = 9, where Tg is the 
temperature in units of 10® K. This temperature is still 
high enough for the system to be the nuclear statistical 
equilibrium (NSE). At that point, the fluid material is 
mainly composed of free nucleons under NSE. Therefore, 
the initial composition of the fluid material is taken to 
be Xp = Ye^g and = 1 — Ye, 9 , where Ye_g is the ini¬ 
tial electron fraction for the nucleosynthesis calculation, 
which is set as the electron fraction when the tempera¬ 
ture of the wind becomes Tg = 9. 

We use a nuclear reaction network for our nucleosyn¬ 
thesis calculations. The network co nsists of 5406 iso¬ 
topes (see Table |3t extended from lYoshida fc Umed^ 
(j20lHl and iSuzukiet al.l 1)201(1( 1. We include all rele¬ 
vant reactions, i.e., thermonuclear reactions and their 
inverse reactions, and weak interactions such as /?-decays 
and electron or positron captures. Thermonuclear reac- 
tion rates are t aken from the JINA REACLIB database 
(iCvburt et al.l 2010(1. W e include the expression of 
iGraboske et al.l ( 1973( 1 and lltoh et all ()1979(1 for the cor¬ 
rection for thermonuclear reaction rates due to the elec- 
tron screening. Weak react i on ra t es are taken from 
Langanke fc Martmez-Pinedol (|2001(1 . lOda et ^ ()1994ll . 


Takahashi fc Yokoil ()1987ll . and FFN85 for the reactions 


whose dependence on the temperat u re an d the density 
is given. We follow iHoriguchi et al.l (|1996ll for weak re- 



Figure 1. Radial profiles of entropy (top), velocity (middle), and 
temperature (bottom) for the two wind trajectories. The wind 
of the mPNS model is shown by solid lines and the wind of the 
ordinary PNS is shown by the dashed lines. The important radii are 
marked with symbols, including the critical point (diamond) and 
the radius of alpha-particle formation (the point of T = 0.5 MeV = 
5.8 X 10® K, circle). 

action rates whose dependence on them is not given, al¬ 
though the experimental rates are known. For the nuclei 
that lack experimental decay-rate information, we use 
Tachiban a (2000, private com munication) for /3+-decay 
rates and iMoller et al.l (12003(1 for /3“-decay rates. In ad¬ 
dition, we include neutrino absorption reactions of free 
nucleons using Equations (1^ and (1^ . 

4. RESULTS 

4.1. The Properties of the Wind Trajeetories 

Eirst, we discuss the properties of constructed winds. 
In Eigure [U we show the contours of the electron frac¬ 
tion at the beginning of the nucleosynthesis calculation 
(Tg = 9) Ye_g, the entropy s at T = 0.5 MeV, and the 
temperature-decrease timescale at T = 0.5 MeV, tt.o.s, 
which are important for the nucleosynthesis. 

The electron fraction at the beginning of the nucle¬ 
osynthesis calculation Ye^g, which is set from the wind 
solutions, has a wide variation from Ri 0.40 — 0.65 in our 
parameter region. It becomes 0.50 when is about 4 
MeV smaller than A larger difference between 
and provides Y^^g < 0.5, i.e., neutron-rich wind, and 
a smaller difference provides proton-rich wind. As shown 
in the next subsection, the electron fraction strongly af¬ 
fects the abundance distribution of the wind material. 
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Table 2 

A List of Nuclei in Our Nuclear Reaction Network. 


Element 

A ■ 

.^min 

A 

.^max 

Element 

A ■ 

.^min 

A 

.^max 

Element 

.^min 

A 

.^max 

Element 

.^min 

.^max 

Element 

A ■ 

.^min 

.^max 

n 

1 

1 

Ar 

27 

67 

Kr 

63 

124 

Xe 

103 

182 

Hf 

149 

240 

H 

1 

3 

K 

29 

70 

Rb 

66 

128 

Cs 

106 

185 

Ta 

151 

243 

He 

3 

6 

Ca 

30 

73 

Sr 

68 

131 

Ba 

108 

189 

w 

154 

247 

Li 

6 

9 

Sc 

32 

76 

Y 

70 

134 

La 

110 

192 

Re 

156 

250 

Be 

7 

12 

Ti 

34 

80 

Zr 

72 

137 

Ce 

113 

195 

Os 

159 

253 

B 

8 

14 

V 

36 

83 

Nb 

74 

140 

Pr 

115 

198 

Ir 

162 

256 

c 

9 

IS 

Cr 

38 

86 

Mo 

77 

144 

Nd 

118 

201 

Pt 

165 

260 

N 

11 

21 

Mn 

40 

89 

Tc 

79 

147 

Pm 

120 

205 

Au 

167 

263 

o 

13 

22 

Fe 

42 

92 

Ru 

81 

150 

Sm 

123 

208 

Hg 

170 

266 

F 

14 

26 

Co 

44 

96 

Rh 

83 

153 

Eu 

125 

211 

Tl 

173 

269 

Ne 

15 

41 

Ni 

46 

99 

Pd 

86 

156 

Gd 

128 

214 

Pb 

175 

273 

Na 

17 

44 

Cu 

48 

102 

Ag 

88 

160 

Tb 

130 

218 

Bi 

178 

276 

Mg 

19 

47 

Zn 

51 

105 

Cd 

90 

163 

Dy 

133 

221 

Po 

182 

276 

A1 

21 

51 

Ga 

53 

108 

In 

92 

166 

Ho 

136 

224 

At 

186 

279 

Si 

22 

54 

Ge 

55 

112 

Sn 

94 

169 

Er 

138 

227 




P 

23 

57 

As 

57 

115 

Sb 

97 

172 

Tm 

141 

230 




S 

24 

60 

Se 

59 

118 

Te 

99 

176 

Yb 

143 

234 




Cl 

26 

63 

Br 

61 

121 

I 

101 

179 

Lu 

146 

237 





Note. 

^He, ®Be, and are not included in the network. 



Figure 2. Contours of the electron fraction, the temperature- 
decrease timescale, the entropy, and the mass outflow rate as func¬ 
tions of Ci/g and . The other parameters are adopted from Table 
^ The black solid lines indicate the electron fraction at the begin¬ 
ning of the nucleosynthesis calculation ^ 6 , 9 - The green dotted lines 
indicate the temperature-decrease timescale, rx- at T = 0.5 MeV in 
units of ms. The red dashed lines indicate the entropy per baryon 
s at T = 0.5 MeV in units of The blue dot-dashed lines indi¬ 
cate the mass outflow rate M in units of 10 “^Mq s“^. The region 
colored in blue is the region where ye ,9 < 0-5, and the region in or¬ 
ange is the region where ye ,9 > 0.5. The regions where the electron 
fraction deviates from 0.5 have deeper color. The shaded region is 
the region where and thus is an unrealistic parameter 

range. 

The obtained electron fraction is close to the value Ve.eq 
from the balance of the two reactions: p+Pe n+e'^ and 
n + i/g p + e~. The value Yr.rq is given appr oximately 
from Equation (77) in lQian fc WooslevI (1199611 . 

The entropy of the wind is distributed in the range 
s « 125 — 145. It gets smaller as the average energy 
of the neutrinos gets larger. This can be understood 
from the equation of the entropy change, i.e., Equation 
(1^ . As the average energy of neutrinos becomes larger, 
the heating rate and T become larger, so the entropy 
increase due to the heating by neutrino absorption is 
suppressed by the high temperature. On the other hand, 
AQ is determined by the gravitational binding energy of 


PNSs and is almost independent of the average energy of 
neutrinos. Thus, the entropy is smaller when the average 
energy of neutrinos is larger. 

The temperature-decrease timescale at T = 0.5 MeV, 
'Tt,o. 5 , which is defined as Equation (1^ . slightly de¬ 
creases (from 6 to 7 ms) as the average energies of the 
neutrinos decrease. This can be understood using the 
last line of Equation (1501) . It is sufficient to pay attention 
to the first term, which is proportional to u/r|T=o.5 MeV- 
This is because the radius where the temperature be¬ 
comes 0.5 MeV is larger than the radius of the critical 
point and thus the second term, u~^du/dT = du/dr, is 
subdominant compared to the first term. The radius 
where the temperature becomes 0.5 MeV gets larger as 
or becomes large because of more efficient heating. 
On the other hand, the terminal velocity is independent 
of or because it is determined by the gravitational 
binding energy at the launching point of the fluid. There¬ 
fore, the winds having small Ep^ or Ep^ have small 
because they have large r|7'=o.5 MeV- 

The mass outflow rate M is on the order of 
10 “^Mq s“^ throughout the figure and becomes large 
as the average energy of neutrinos is large because of the 
efficient heating. 

The contours of the temperature-decrease timescale, 
the mass outflow rate, and the entropy are almost par¬ 
allel to each other. This is because all these quantities 
depend mainly on the heating rates, i.e., on Ep^ and Ep^ 
(see Equations (fHl) . (fT71) and (flSll l. On the other hand, 
the electron fraction Ye ,9 depends on the difference of 
the average energies between electron neutrino and elec¬ 
tron antineutrino — Ep_ (from the appr oximate form 
of Equation (77) in lQian fc Woosle'^ll996[) . 

The average energies of electron neutrinos and electron 
antineutrinos usually have the relation Ep^ < Ep^, and the 
shaded region in Eigure[2]is physically less relevant. This 
is because PNSs have much neutrons than protons, so 
the radius of the neutrinosphere of electron neutrinos is 
larger than that of electron antineutrinos. 

4.2. Dependence of the Abundance Distribution on the 
Average Energies of Electron Neutrinos and 
Electron Antineutrinos 
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Figure 3. The contours of the average mass number of final 
abundance distribution (A), the electron fraction at the beginning 
of the nucleosynthesis calculation Ye, 9 , and the entropy s. The 
solid lines indicate the average mass number (A) of the heavy- 
element set. The red doted lines indicate the electron fraction 
Ye, 9 . The green dashed lines indicate the entropy s per baryon at 
T = 0.5 MeV in units of k^. The r-process and the i^p-process 
proceed in blue and orange region, respectively. The deeper the 
color of a region is, the larger mass number nuclei are produced in 
the trajectory in the region. The shaded region at the lower right 
of the black-dashed line covers unrealistic parameter region where 

Ci/g > Cj/g . 


We calculate the nucleosynthesis for the trajectories 
with parameters listed in Table [H In Figure [3l three 
contours for the electron fraction at the beginning of 
the nucleosynthesis calculation Y^^g, the entropy s at 
T = 0.5 MeV, and the average mass number of the syn¬ 
thesized heavy elements {A), defined below, are shown 
in the ~ plane. In order to see the typical mass 
number of heavy elements synthesized in each trajectory, 
we define the heavy-element set M >48 as the nuclei that 
have the mass number A > 48, which includes iron-group 
and heavier elements, and the average mass number of 
the heavy-element set of the final abundance distribution 
{^) as 


A) = y^ E a,y{z,,a,) 


(31) 


^€v4>48 


where Y{Zi,Ai) is the abundance in number of the 
species having the atomic number Zi and the mass num¬ 
ber Ai. In Equation (ISTI) . Yh is the abundance of the 
heavy-element set defined as 

Yh= Y. (32) 

A>48 

Figure [3] is divided into two regions parted on the con¬ 
tour of Ye _9 = 0.5. In the region where Y^^g < 0.5 
(the blue region in Figure [3]), heavy elements are syn¬ 
thesized through the r-process, while in the region where 
Ye ,9 >0.5 (the orange region in Figure [3|), the t'p-process 
is the main process that produces heavy elements. 

The average mass number (A) strongly depends on 
the electron fraction. In the region where |Ye ,9 — 0.5| 
is larger, heavier elements are synthesized. In the region 


with Ye ,9 < 0.45, the average mass number exceeds 130 
through the r-process. On the other hand, strong vp- 
process proceeds in the region with Ye ,9 > 0.55 and the 
average mass number becomes (kl) > 80 — 100. For the 
synthesis of heavy elements, higher entropy is favored 
because high entropy prevents seed elements from being 
produced by strong photo-disintegration and the number 
of neutrons and protons captured by each seed nucleus in 
the subsequent phase becomes large. We see that even 
for the same electron fraction, the trajectory that has 
higher entropy has larger average mass number for the 
above reason. However, the dependence of the average 
mass number on the entropy is weaker than the depen¬ 
dence on the electron fraction. 


4.3. Details of the Nucleosynthesis Processes 

In the previous subsection, we showed that heavy el¬ 
ements with A > 100 can be synthesized through the 
r-process or the zzp-process in the neutrino-driven winds 
of the mPNS. Here we show more details for the aver¬ 
age mass number of the heavy-element set and the main 
products in nine wind trajectories (a)-(f) in Figure [S] In 
Table m we list the properties of these wind trajectories 
and the nucleosynthesis results of them. Trajectory (a) 
is the wind with the fiducial parameters (see Table [l]). 
In order to see the effect of the entropy, we consider tra¬ 
jectories (b) and (c), which have almost the same Ye ,9 as 
trajectory (a) but different s (130— 138). Trajectory (d) 
has almost the same s as trajectory (a) but lower Yg 9 of 
0.42. 

In order to see the nucleosynthesis in the proton-rich 
region, we choose trajectory (e), which has the same en¬ 
tropy as trajectory (a) but a higher electron fraction of 
0.55. Trajectories (f) and (g) have the same Ye ,9 as tra¬ 
jectory (e) but different s (129 — 141). We choose tra¬ 
jectory (h) to see the dependence of Ye, 9 . As an extreme 
case, we choose trajectory (i), which has the largest (A) 
in the physically relevant proton-rich region. 

In order to see gross features of the heavy-element syn¬ 
thesis, we define the average rates of individual reactions 
per nucleus in the heavy-element set >l> 48 , defined in the 
previous subsection, as 




Yu 


E pNA{a^p,-,){Z,,A,)v)YpY{Z^,A,), (33) 

i£A>48 




Yh 


E 

iGA>48 


^{■y,p){Zi -\-l,Ai-\- l)Y{Zi -|- 1, Ai -|-1), 


(34) 




Yh 


E P^A{(^(n,j){Zi,A,)v)YnY{Zi,Ai), 

^GA>48 


(35) 


A(7,n)— ^ E ^p,n)[Zi,Ai-\-\)Y{Zi,Ai-\-l), (36) 

^ i^A>A8 

\n,p) = ^ E P^A{c^in,p){Zi,A,)v)YnY{Zi,A,), (37) 


\p,’n) = Y E P^A{cr(p,n){Zi-l,Ai)v)YpY{Zi-l,Ai), 

^ ^GA>48 

(38) 











nucleosynthesis in neutrino-driven wind in HNe 


9 


Table 3 

Properties of Trajectories 



(MeV) 

Ve,9 

'^T,0.5 

(ms) 

s 

(10 

M 

-^Mqs 

0 

0 

^71 

(T) 

Mainly Produced Nuclides 

(a) 

11.0 

16.0 

0.481 

6.59 

133 

4.62 

0.696 


118 


(b) 

9.29 

14.2 

0.482 

6.35 

138 

3.66 

0.730 


123 

125,126,128,130rpg 127j 129,131,132,134 

(c) 

12.5 

17.4 

0.480 

6.76 

130 

5.54 

0.676 


115 

109Ag,lll'116Cd, 122,124g„^ 121.123gb^ 125,126^6, 127i 

(d) 

9.29 

17.7 

0.421 

6.60 

133 

4.67 

0.795 


139 

134Xe, 133Cs, 135.137Ba^ 177,178,179Hf, 198pt, 198 Au 

(e) 

13.1 

14.2 

0.545 

6.59 

133 

4.63 


39.7 

92.4 

96,98,100j:j^y 102,104p^ 10®Cd 

(f) 

10.2 

12.2 

0.543 

6.23 

141 

3.33 


49.2 

101 

98,100j^^ 102,110pj^ 1*^9 Ag 108,lll,113Q(i 

(g) 

15.4 

16.0 

0.545 

6.84 

129 

5.95 


31.1 

84.4 

84g,,^ 102,104p^^ 108(g^ 

(h) 

13.7 

13.6 

0.566 

6.59 

133 

4.60 


55.5 

106 

98,100j:^y llOpd 108,111,113,116q^^ 

(i) 

11.0 

11.0 

0.587 

6.18 

141 

3.21 


110 

129 

isSLa 


Note. In the column “Mainly Produced Nuclides,” we list top 10 nuclides which have the production factors larger than 1/10 times the 
largest production factor in each trajec tory . The nuclide which has the largest production factor in each trajectory is written in bold. The 
parameter /200 is defined as Equation ll43tl . 


^ iGA>48 


Xp- - 


1 

Yh 


Y, Xp-{Z,,A,)Y{Z,,A,). 

*£^>48 


Here ^ Ji.'j^ , nnd 

<^{n,p){Z^ A) are the cross sections of (p, 7 ), 

{p,n), and {n,p) reactions of nuclear species "^Z. 
X^y p'){Z, A) and A) are reaction rates of the 

photo-disintegrations ( 7 ,p) and ( 7 ,n). Xp+{Z,A) and 
Xp-{Z,A) are /3+ decay and /3“ decay rates. The rates 
defined in Equations dssD-dini) represent the reaction 
rates of dominant isotopes at a given time. In addition, 
we define the average proton number {Z) of the heavy- 
element set and the average neutron number (N) of 
them as 


(Z) 

(N) 


Yh 


Yh 


Y ZAY(,Z,,A,\ 


^€v4>48 


^ N,Y{Z,,Ai). 


4 €- 4>48 


(41) 

(42) 


4.3.1. Nucleosynthesis Result for the Fiducial Parameter Set 

(a) 

First, we show the result of the nucleosynthesis of the 
fiducial wind trajectory with = 11 MeV and = 
16 MeV, the parameter s et of which is based on the result 
of ISekiguchi et al.l (|2012f) . This trajectory corresponds to 
point (a) in Figure O The electron fraction of this 
trajectory is 0.48. The r-process occurs, and the average 
mass number becomes (A) = 118 (see Table [3|. 

Figure 0] shows the time evolution of the mass frac¬ 
tions of neutrons, protons, and the heavy-element set 
(top panel) and the average reaction rates of the heavy- 
element set (A^p^^p A^.y pp 

Xp-, and Xp+; bottom panel). We show the snapshots of 
the abundance distribution in Figure [SJ 

We explain the time evolution of the nucleosynthesis 
in this trajectory. Nucleons are still free particles at a 
temperature higher than Tg « 6 . After the temperature 
drops below Tg w 6 , the heavy elements up to A ^ 100 
are temporarily produced, and the average charge num¬ 
ber becomes (Z) = 38 (see top panel of Figure [Hand [5]- 
(1)). This is because the numbers of nucleons are larger 
than that of the corresponding quasi-statistical equilib- 



0.001 0.01 0.1 1 
time[s] 


Figure 4. Top: time evolution of the mass fractions of protons, 
neutrons, and the total value of the heavy elements, and the aver¬ 
age proton and neutron num ber of them in trajector y (a), which has 
the parameter set based on ISekiguchi et al.l 1120121 V Bottom: time 
evolution of the average reaction rates per nucleus in the heavy- 
element set in the same trajectory. The rates of ( 11 , 7 ), (p, 7 ), and 
(n,p) reactions are shown as green, red, and blue solid lines, re¬ 
spectively, and their inverse reaction rates are shown in dashed 
lines. The rates of /3~-decays (solid line) and /3'*^-decays (dashed 
line) are shown as pink lines. In addition to them, the temperature 
is shown as a black solid line, and the temperature-decrease rate 
Ty ^ is shown as the black dashed line. 


rium (QSE) value ()MeveHl2002[ ). After this phase, the nu¬ 
clei heavier than iron-group are photo-disintegrated into 
lighter nuclei as the free nucleons are recombined into 
alpha particles at Tg ri 5.5. The nuclei of A ^ 80 — 90 
are not disintegrated completely in this phase because 
the nuclei at the neutron magic number at V = 50 have 
relatively larger binding energies (see Figure [S]- ( 2 )). 

The r-process begins after the temperature decreases 
below Tg « 4. In this phase, neutron capture reactions 
proceed until heavy nuclei reach the classical waiting 
points, where the reaction (n, 7 ) balances with ( 7 , n) (see 
top panel of Figured and Figure [5} (3)). In the ordinary 
r-process, the nuclei on the waiting points can obtain 
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additional neutrons through the /3“ decay. In the case 
here, however, the {p, n) reactions play the role. The 
protons are supplied through the electron neutrino ab¬ 
sorption reaction of free neutrons [n + Ve —>■ p + e~). 



t = 2.88x10 s 
Tg= 3.690 
p = 1.21 xIO^ g cm'" 


1 = 7.69x10'^s 
T,= 0.497 
p = 1.35xIO gem'" 

150 



-14 


200 


Figure 5. Snapshots of the nuclear abundance distribution of 
nucleosynthesis in the N — Z plane for fiducial trajectory (a). The 
blue dashed lines indicate the nuclei of A = 56, 130, and 195. 
Stable nuclei are shown by black squares. Vertical and horizontal 
belts represent the N = S2 neutron magic number and the Z = 50 
proton magic number, respectively. 


At the time t ^ 0.01 s, the rate of (p, n) reactions 
decreases because of the low temperature and the heavy- 
element production is temporarily frozen at the N = 82 
magic nuclei. After this phase, as the temperature de¬ 
creases further, the rate of ( 7 , n) reactions becomes too 
small to balance with (n, 7 ), so the neutron capture reac¬ 
tions produce the heavy elements beyond TV = 82. In this 
phase, the so-called cold r-process is realized, in which 
(n, 7 ) reactions and /3“ decay react ions have mai n con¬ 
tributions to produce heavy nuclei (jWanaioll2007D . The 
nuclear abundance distribution, as functions of neutron 
and proton numbers, at Tg ss 0.5 is shown in Figure[S]-(4). 
In this figure, we see that, although the fraction is very 
small, nuclei of V > 82 are produced by the r-process. 
After the phase (Tg < 0.5), all reaction rates get small 


and the abundance distribution is completely frozen. 

The middle panel of Figure | 6 ] shows the abundance 
distribution of trajectory (a). The nuclei in the second 
peak elements are mainly produced. The mass number 
of the most abundant nucleus in the heavy element is 
A = 126. 

The middle panel of Figure [7] shows the production 
factors of final abundances of the nucleosynthesis re¬ 
sult of trajectory (a). The production factor is defined 
as the abundance ratio to the solar system composition 
(|Anders fc Greves^ll989H . It is seen that the nuclei that 
form the second-peak are mainly produced (see mainly 
produced nuclides in Table E]). 

4.3.2. Parameter Sets in Neutron-rieh Region {b)-{d) 

Next, we show the dependence of the abundance dis¬ 
tribution on the neutrino parameters in the trajectories 
in the neutron-rich region Ye^g < 0.5 ((b)-(d)), where the 
r-process occurs as in the fiducial trajectory (a). 

The top and bottom panels of Figures El and [7] show 
the abundance distributions and the production factors 
of trajectories (b) and (c), which have the same electron 
fractions as the fiducial trajectory (a) (0.48) and different 
entropies (138 and 130 for trajectory (b) and (c), respec¬ 
tively). We see that the peak of the abundance distribu¬ 
tion of trajectory (b) is A = 129, which is slightly larger 
than that of the fiducial trajectory (a) {A = 126), and 
the second-peak elements such as Te, I, Xe, and Cs are 
produced (see Table [3] for individual isotopes). 

In trajectory (c), on the other hand, the abundance 
peak is A = 125, which is almost the same as trajectory 
(a), and mainly produced elements are also the same (i.e., 
the isotopes of Sn, Sb, Te, and I). In addition to these 
elements, lighter elements such as Ag, Cd, and Pd are co¬ 
produced. This is because the r-process does not proceed 
sufficiently and the nuclei that do not pass through the 
N = 82 magic number finally decay into these elements. 

The abundance distribution of trajectory (b) is slightly 
shifted to heavier-element production, and the distribu¬ 
tion of trajectory (c) indicates an opposite trend, but 
their difference is small. 

We will show the nucleosynthesis result of a more 
neutron-rich trajectory with the entropy similar to tra¬ 
jectory (a). The electron fraction of trajectory (d) is 
Ye,g = 0.42 (see Table[3]). Figure [ 8 ] shows the abundance 
distribution and the production factors of the trajectory. 
We clearly see the second peak and third peak in the dis¬ 
tribution. This means that a stronger r-process occurs 
and the nuclei of wide mass number range from second 
peak to the third peak are produced. The time evolution 
of the average reaction rates is almost the same as in tra¬ 
jectory (a) but the neutron-to-seed ratio at the beginning 
of the cold r-process is about 30. Therefore, we conclude 
that the nucleosynthesis results do not depend strongly 
on the entropy, but mainly on the electron fraction Ye.g- 

The combination s/Ye^gr^^^ ^ is a useful quan¬ 
tity for estimating the neutron-to-seed ratio and 
whether the r-proces s produces the third-peak elements. 
iHoffman et al.l ()1997t) gave an approximate condition for 
making the third-peak (A ^ 200) elements in the trajec¬ 
tory as 

^ (^/230) > 

" (Ye. 9 / 0 . 4 )(TT, 0 . 5/20 ms)i/3 ) 
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Figure 6. Abundance distributions of trajectories (b) (top), (a) 
(middle), and (c) (bottom). These trajectories have almost the 
same electron fractions and different entropies (see Table 0. 



Figure 8. Final abundance distribution (top) and the production 
factor (bottom) of trajectory (d). 

(This definiti on of /200 is based on Equation (1) in 
IWanaicill2013L .) Note that the condition is somewhat 
relaxed for the wind where Ye > 0.46 (see Figure 10 in 
I Hoffman et al.l 119971 ) because the rate of seed element 
production decreases in such a high-Yg circumstance. 

The values of /200 are shown in Table [31 It is about 
0.696 in trajectory (a), which is below unity. As a re¬ 
sult, the neutron-to-seed ratio is not large, 10-20 at the 
beginning of the cold r-process, so the r-process in this 
trajectory is not strong enough to produce the third-peak 
nuclei. On the other hand, the values of /200 are 0.730 
and 0.676 for trajectories (b) and (c), respectively, so 
the neutron-to-seed ratio is bigger in trajectory (b) and 
smaller in trajectory (c). Therefore, slightly heavier and 


Figure 7. Production factors of the nuclei obtained from the 
nucleosynthesis in trajectory (b) (top), (a) (middle), and (c) (bot¬ 
tom). The solid line connects isotopes of a given element. The 
dotted horizontal lines indicate the largest production factor and 
a factor of 10 smaller than that. The dashed horizontal line indi¬ 
cates the median value. 

less heavier nuclei are synthesized in trajectories (b) and 
(c), respectively, as shown in Figure |6l 

In trajectory (d), /200 = 0.795, which is larger than 
those of the above trajectories because of the smaller 
value of Yg. Therefore, a stronger r-process proceeds and, 
as shown in Figure |S1 the third-peak elements are pro¬ 
duced in the trajectory. (The third peak appears while 
the value /200 is still smaller than the unity. This is be¬ 
cause Equation (l43)l is only an approximate condition.) 

4.3.3. Parameter Sets in Proton-rich Region (e)-(i) 

The parameters suggested in iSekieruchi et al.l ()2012[ ) 
are close to the contour of Yg.g = 0.5, so the initial elec¬ 
tron fraction Yg^g of the trajectories in a certain range of 
the parameter values and becomes greater than 
0.5. In these trajectories, the r'p-process occurs. Here we 
explore the nucleosynthesis and the abundance distribu¬ 
tions of some proton-rich trajectories (e)-(i). 

Figure m shows the time evolution of the nucleosynthe¬ 
sis and the abundance distribution of trajectory (e), of 
which Yg 9 and s are 0.55 and 133, respectively. Initially, 
the material is composed mainly of free protons and neu¬ 
trons. Then, heavy nuclei are temporarily produced as 
in the neutron-rich case (see top panel of Figures |9| and 

mil)). 

At temperature Tg ss 4.3, the heavy nuclei are com¬ 
pletely photo-disintegrated into iron-group nuclei and the 
mass fraction of neutron A„ decreases to ^ I 0 “® (see 
top panel of Figures [9] and [T0l(2)). After this phase, 
neutrons are supplied through the electron antineutrino 
absorption on free protons (p -b T'g —>■ n -I- e+). After the 
temperature decreases to Tg Ri 4, the balance between 
(n,p) and {p,n) reactions is broken, and then the vp- 
process starts to operate. In this phase, (p, 7 ) and ( 7 ,p) 
reactions are balancing, and most of the heavy nuclei are 
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Figure 9. Same figure as Figure[4l but for the result of trajectory 

(e) . 

in the waiting point of the ^'p-process. The production of 
the heavy elements occurs via (n,p) reactions. The path 
of the i^p-process ends at Z = 50 proton magic nuclei. 
The reaction flow hardly passes through Z = 50 because 
of small cross sections of (p, 7 ) reactions of Z = 50 iso¬ 
topes (see Figure [101-1311. When the temperature gets 
below Tg Ri 2.5, the balance between ( 7 ,p) and (p, 7 ) 
reactions in broken. In this phase, (n,p) reactions grad¬ 
ually take over from ( 7 ,p) reactions since neutrons do 
not feel the Coulomb barrier. 

When the temperature decreases below Tg ss 1.5, the 
heavy nuclei are concentrated at the stable line and (n, 7 ) 
reaction rates become larger than (n,p) reaction rates. 
As a result, a fraction of large A nuclei become neutron- 
rich (see Figure fTUl-fd')'). After this phase, all reaction 
rates get small and the abundances are completely frozen. 

The abundance distribution of trajectory (e) (middle 
panel of Figure [TT]) indicates a broad peak of A ^ 60 — 
120. The middle panel of Figure IT^ shows the production 
factors of final abundances of this trajectory. It is seen 
that A ^ 100 proton-rich isotopes, such as Ru, Pd, and 
Cd, are produced. 

We also consider trajectories (f) and (g), of which the 
electron fraction is almost the same as for trajectory (e) 
but the entropy is different (141 and 129 for trajectories 

(f) and (g), respectively). The top and bottom panels 
of Figure [TT] show the abundance distributions of these 
trajectories. The r'p-process scarcely passes through the 
Z = 50 magic number in either of the trajectories, so 
the peaks of their abundance distributions are almost 
unchanged at A 100. The top and bottom panels of 
Figure [T2| show the production factors of these trajecto¬ 
ries. They are almost the same as that of trajectory (e) 
except that the peak of the production factors slightly 
shifts to a larger mass number in trajectory (f). 

In trajectories that have larger Ve^g, heavier nuclei are 
produced compared with trajectory (e). In trajectory 



Figure 10. Snapshots of the nuclear abundance distribution of 
nucleosynthesis for the trajectory (e). 

(h), which has the electron fraction Ye ,9 = 0.57, for ex¬ 
ample, the isotopes of Ru, Pd, Cd, In, and Sn are abun¬ 
dantly produced (see Table [3] for details). However, the 
produced elements that have A > 100 are not proton- 
rich because of efhcient (n, 7 ) reactions in the phase of 
Tg<2. 

In trajectory (i), the most extreme (but not unphysi¬ 
cal) trajectory in FigurejS] efhcient (p, 7 ) reactions due to 
a large number of free protons proceed strong i^p-process 
nucleosynthesis and heavy nuclei can be produced be¬ 
yond the Z = 50 magic number. Figure [T3| shows the 
nuclear abundance distribution at Tg w 0.55 in trajec¬ 
tory (i). It is seen that a large amount of nuclei pass 
through Z = 50 and are concentrated in a neutron-rich 
region. 

In Table 121 we show the quantity A„, the number of 
neutrons per seed nucleus produced via p + Vp, ^ n + e~*~ 
react ion for Tg < 3 (see Equation (1) in iWanaio et al.l 
Unni). When A„ is larger, more (n,p) reactions occur 
during the t'p-process so that heavier nuclei are pro¬ 
duced. In higher-entropy trajectories, less seed nuclei are 
produced and A„ is larger. In trajectory (i), A„ = 110, 
which is much larger than that of (e)-(g) (39.7, 49.2, 
and 31.1, respectively). This is the reason why the vp- 
process can pass through the Z = 50 magic number and 






































































nucleosynthesis in neutrino-driven wind in HNe 


13 




mass number A 


Figure 11. Same as Figure^ but for the proton-rich trajectories 
(f) (top), (e) (middle), and (g) (bottom). 



Figure 13. Nuclear abundance distribution at Tg ^ 0.55 in tra- 
jectory (i). 



Figure 14. Final abundance distribution (top) and production 
factor (bottom) of trajectory (i). 

the nuclei are concentrated in a neutron-rich region in 
t his trajectory. 

IWanaio et al.l (|201ir) , a comprehensive study of the vp- 


Figure 12. Same as Figure[7l but tor the proton-rich trajectories 
(f) (top), (e) (middle), and (g) (bottom). 

process in the ordinary PNS wind, showed that heavy 
p-nuclei of A > 110 can be synthesized through the vp- 
process only in extreme conditions where Yg > 0.6. In 
g- our HN wind model, however, the i^p-process produces 
^ heavy nuclei of ^ > 100 for smaller values of Yg, although 
ra these are not p-nuclei. This is because the higher neu¬ 
trino luminosity causes less seed element production and 
more efficient neutron production through p-l-i'e n-l-e"'" 

reactions. As a result, the value of A„ can become larger 
than 100 , which is not ach ieved even if Yg is s et as 0.8 
in the ordinary PNS wind (jWanaio et al.l 1^1 Ih . Heavy 
p-nuclei of A = 110 — 125 ar e underproduced in the 7 - 
process l|Rauscher et al.l[2002D . so our HN wind may be a 
possible site for this p-nuclei synthesis if produced nuclei 
are p-nuclei. This may be possible when we consider the 
wind termination shock as indicated in Section [S] 

Only a small amount of nuclei in the neutron-rich re¬ 
gion go beyond the TV = 82 magic number during the 
nucleosynthesis. Thus, the A ^ 130 — 140 peak is pro¬ 
duced in the final abundance distribution, as shown in 
the top panel of Figure fT4l The bottom panel of Figure 
[14] shows the production factors of this trajectory. The 
figure shows that ^^®La is mainly produced. 

5. DISCUSSION 

5.1. The Effect of Heavy-element Synthesis on the 
Galactic Chemical Evolution 

In previous sections, we investigated the heavy-element 
synthesis in the neutrino-driven wind from mPNSs in 
the HN scenario. It is suggested that HNe may have 
contributed to the synthesis of heavy elements and to the 
Galactic chemical evolution in the early universe because 
HNe may have occurred in Population HI ('Pop FI) st ars 
(e.g.. lUmeda fc Nomotoll200,5l iTominaga et alJl^(I7alfO ). 
Therefore, here we discuss the effect of heavy-element 
synthesis in the mPNS wind on the Galactic chemical 
evolution. 

It is known, from observations, that a large amount of 
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®®Ni is produced in HNe. Therefore, we discuss the effect 
of HNe on the Galactic chemical evolution by comparing 
the production factors of heavy elements produced in the 
neutrino-driven wind with that of ®®Ni (Fe) produced 
in HNe. If the production factor of a heavy element is 
larger than or equal to that of Fe, it can be said that the 
heavy-element synthesis in the neutrino-driven wind of 
HNe affects the Galactic chemical evolution. 

We will evaluate the total production factors of ele¬ 
ments considering the total ejecta of HNe. The total 
production factor of Fe can be evaluated from observa¬ 
tional analyses of the ®®Ni yield in recent HNe. The total 
production factor of Fe, Ptot(Fe), is evaluated as 


-Ptot(Fe) 


M(56Ni) 

YQ(Fe)Mi„i 


(44) 


where M(®®Ni) is the ®®Ni yield, Y 0 (Fe) is the solar mass 
fraction of Fe, which is « 1 x 10“^, and Mini is the initial 
mass of the HN progenitor. TablelHshows the properties 
of observed HNe and the total production factor of Fe. 
We see from the table that the total production factor 
of Fe Ptot(Fe) is about > 4. Thus, the total production 
factor of the element that can contribute the Galactic 
chemical evolution should be greater than ^ 10 at least. 
The total production factors of the heavy elements pro¬ 
duced in the wind are estimated as 


P,n,(^Z) = ^^^l^Pwind(^^) 


MwindTNS 


M, 

G0-^Pwind(^^) 


-^ej ,tot 
^-4 td (A > 


Pwmd(^^) 


Mwind 


4 X 10~^Mq S“1 


-1 


TNS I ( Mejyot j 

T7 j 1^4OM0 j ’ 


(45) 


where Pwind("^^) is the production factor of the element 
in the wind ejecta, which is shown in Section|Tl M^ind 
is the mass outflow rate of the wind, and tns is the life¬ 
time of mPNSs. Thus, the elements that have the pro¬ 
duction factor of winds Pwind("^-^) ^ 10 ^ may contribute 
to the Galactic chemical evolution. 

For the neutron-rich wind near the fiducial parameter 
set, as seen in Figure [71 the second-peak elements such 
as Sb, Sn, Te, and I are co-produced with these produc¬ 
tion factors larger than about 10®. Therefore, the mPNS 
wind may contribute to the Galactic evolution of these 
elements. In a more neutron-rich region, e.g., in trajec¬ 
tory (d), the wide range of nuclei from the second peak 
to the third peak have sufficient production factors to 
contribute to the Galactic chemical evolution as seen in 
Figure [TH 

On the other hand, in the proton-rich wind, proton- 
rich isotopes of H ^ 100, such as Ru, Pd, and Gd have 
sufficient production factors larger than about 10 ® (see 
Figure [12] for trajectories (e)-(g)). In more proton-rich 
wind, e.g., trajectory (i), the production factor of ^®®La 
is ~ 10 ^, which is much larger compared with other ele¬ 
ments produced in this trajectory as seen in Figure [Td] 

We see from Table [4| that the total production factor 


of Fe tends to be larger when the progenitor mass of the 
HN is larger (for SN 1999as, it is about 50 at least). 
Thus, the threshold of the wind production factor for 
contribution to the Galactic chemical evolution may be 
one order larger than ~ 10® for the HNe of more massive 
stars. 

5.2. Heavy-element Abundance in Very Metal-poor Stars 

Gan we explain the abundance patterns of r-rich, very 
metal-poor stars by this mPNS wind model? As seen 
in the previous sections, the so-called main r-process, 
which produces all r-process elements up to the third 
peak, hardly occurs in our parameter sets. Therefore, it 
is very hard to explain the origin of very me tal-poor stars 
having the “universality” of the r-process (jSneden et al.l 
1200.11 [20?)^ using our mPNS wind model. 

There are, however, some very metal poor stars that 
do not have that universality. For example, HD 122563 
is one of the stars, whose abundance pattern shows a 
grad ually decreasing tr end as a function of atomic num¬ 
ber ([Honda et al.ll200'^ . They are called “weak r stars.” 
Thus, we compare our result with the abundance pattern 
of these weak r stars. Figure [15] shows the abundance 
pattern of the nucleosynthesis result of the trajectory just 
near point ( a) in Figure [3] which h as the parameter set 
indicated in iSekignchi et al.l (|2012t) . The parameter set 
of this trajectory is = 11.0 MeV, = 16.8 MeV. 
In this figur e, the scaled obse rvational abundances of 
HD 122563 ([Honda et al.l 120(3^ are shown with points. 
The abundance pattern shows an approximate agreement 
with the observational trend in Z = 39 — 45 and Z > 58. 
Recent observations show that the heavy-element abun- 
dance pattern of weak r stars has som e varieties (e.g., 
iHonda et al.ll200/i: iRoederer et al.l[MT(ill . These varieties 
may be explained by considering the difference of the 
properties of mPNS winds. 



Figure 15. Abundance distribution of the elements as a function 
of atomic number Z. The solid line indicates the nucleosynthesis re¬ 
sult of the trajectory which has the parameter set Ci/g = 11.0 MeV, 
= 16.8 MeV. The das hed line is the scaled solar abundances 
JAnders fc Grevessd 1198911 . Th e points with error bars are the 
scaled abundances of HD 122563 dHonda et al.l2006ll . The inverted 
triangles show the upper limits of the elements. The observational 
abundance patterns are scaled to the Eu abundance of the nucle¬ 
osynthesis result. 


5.3. The Dependence of Heavy-element Synthesis on the 
Mass, Radius, and Neutrino Luminosities of PNSs 
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Table 4 

Properties of Observed Hypernovae 


Name 

Type 

E 51 

44im/440 

M(»»Ni)/M0 

Ptct(Fe) 

Reference 

SN 1997ef 

ic 

8 

30 - 35 

0.15 

4-5 

iwamoto et al. (20UU) 

SN 1998bw 

Ic 

30 

40 

0.7 

17.5 

Iwamoto et al. 11998) 

SN 2002ap 

Ic 

7 

21 

0.10 

4 

Mazzali et ai. (2002) 

SN 2003dh 

Ic 

30 - 50 

35-40 

0.35 

8-9 

Mazzali et ai. (2003) 

SN 1999as 

Ic 

20 - 50 

60 - 80 

> 4 

>42-60 

Denpj et al. (2001) 

SN 20031w 

Ic 

60 

40 - 50 

0.55 

11 - 14 

J)3azSl^t al. (2006) 


Note. f? 5 i is the explosion energy of SN in units of 10^^ erg. 


In this work, we investigated the dependence of heavy- 
element synthesis on the average energies of electron neu¬ 
trinos and electron antineutrinos. Here we discuss the 
dependence of the yields on the other wind parameters. 

In addition to the uncertainties of average energies 
of electron neutrinos and electron antineutrinos , 

there will be uncertainties of the average energy of other 
flavors of neutrinos originated in approximated treat¬ 
ment of neutrino physics. The neutral-current neutrino 
reactions with ^He will affect heavy-element synthesis. 
However, this effect is less important than the charged- 
current neutrino reactions with n and p. Thus, this un¬ 
certainty will have only a minor impact on the heavy- 
element synthesis. Luminosities of individual flavors of 
neutrinos Lp,,, Lp^, and also have uncertainties for 
the same reason. These uncertainties are expected to af¬ 
fect the heavy-element synthesis through the change of 
the wind properties and the neutrino reactions. 

On the other hand, the uncertainty of the radius of 
(the neutrinosphere of) the mPNS i?p originates not only 
in the matter effect on the neutrino opacity but in the 
fact that the EOS of high-density nuclear matter, or the 
mass-radius relation of the NS, is still unknown. 

The mass of a mPNS M is another parameter that de¬ 
pends on the properties of its progenitor, such as initial 


mass, rotation, metallicit y, and so on. A recent study 
of Pop in star formation ( Susal[2M^ iHirano et al.ll2od 


suggested that Pop HI stars of the mass oyer IOOMq can 
be formed. In the final stage of this eyolution, these mas¬ 
sive stars may have a massive iron core more than ZMq. 
Recent observations suggest that the maximum mass of 
a neutron star is more than 2Mq (iDemorest et al.ll201d : 
lAntoniadis et al.l[2013H . so the EOS of the nuclear matter 
needs to be sufficiently stiff in order to support that mass. 
In addition, Pop M stars cannot lose their angular mo¬ 
mentum because of the absence of the mass loss. There¬ 
fore, when these stars collapse, rapid rotation of these 
cores prevents them from collapsing to black holes, in 
addition to their thermal pressures. Thus, the existence 
of mPNSs is quite possible when Pop HI stars collapse. 
In the studies about the nucleosynthes is of Pop III SNe 
and HNe le.g.. lIJmeda fc Nomotoi[200^ . the “mass cut,” 
the remnant mass after the explosion, was set to explain 
the abundance patterns of the nuclei up to iron-peak el¬ 
ements of extremely metal-poor stars. They changed the 
mass cut for different progenitor masses and explosion 
energies and set as « ?>Mq in some cases. Therefore, 
these studies imply that mPNSs may indeed be formed 
during Pop IH HNe and the masses of them may change 
for different progenitor masses and explosion energies. 
To discuss the history of heavy-element synthesis of HNe, 
further study is needed that covers a wide parameter 
range of the mass of HN progenitors. 


Uncertainties and varieties of these parameters can af¬ 
fect the entropy s, the temperature-decrease timescale 
tt, and the electron fraction Ye of the neutrino-driven 
wind and thus the nucleosynthesis processes. In order to 
investigate the yields of the heavy elements of HNe, the 
systematic search of the nucleosynthesis changing these 
parameters is needed. This also will appear in our next 
work. 

In particular, we find an interesting nucleosynthesis 
feature in the parameter space. There is a region where 
{A) is large but the initial electron fraction Yg.g is only a 
little larger than 0.5, i.e., the region around k, 9 MeV 
and ss 13 MeV. In the trajectories in this region, the 
abundances freeze out without the reaction equilibrium 
between the nucleons and alpha-particle. As a result, 
after the temporal production of the A > 100 nuclei, 
photo-disintegration does not proceed completely, and 
nuclei around the N = 50 magic number remain. This 
process can occur in a realistic parameter range, but de¬ 
tailed analysis of the process in these trajectories is be¬ 
yond the scope of this paper. We will make a detailed 
analysis in our next work. 

5.4. Effect of Wind Termination 

As in the case of ordinary SNe, it could happen that the 
mPNS wind terminates at some radius because of the ex¬ 
istence of the stellar material. The wind termination can 
affect the nucleosynthesis because it chang es the thermal 
histo r ies of the fluid elements in the wind (iKuroc a et al.l 
l2008HArcones fc .lankall^fiTTI : iWanaio et ai.ll20111 . Here 
we discuss the effect of the wind termination shock on 
the nucleosynthesis in the previous sections. 

We assume that a termination shock appears at the 
radius i?sh- We apply the relativistic Rankine-Hugoniot 
jump conditions there. The jump conditions are 

Pw'Uw = PsMs, (46) 

= p^hsUg -t- Pg, (47) 

Pw^w^w\/l 4“ - Pg/lgrtg-y^l -f Wg , (48) 

where h = 1+e+P/p is the specific enthalpy, and the sub¬ 
scripts “w” and “s” indicate the quantity of the upper 
stream (wind) and downstream (shocked) material, re¬ 
spectively. We assume that the velocity of the shock po¬ 
sition is negligible compared to the wind velocity ~ 0.3c. 
The details of solving these equations are described in 
Appendix [B] Solving the jump conditions (|46|) - (|4^ . we 
can evaluate the downstream quantities Us and Tg from 
the upstream quantities Mw and Tw (the electron fraction 
is assumed to be unchanged across the shock). Then u, 
T, and Ye profiles in the shocked region (r > i?gh) are 
obtained by Equations O-® using the shocked values 
Us and Tg as the boundary condition at r = Rgh- 
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5.4.1. The Effect on the r-process 

Here we discuss the effect of the wind termination 
shock on the nucleosynthesis result in the neutron-rich 
region in the parameter space. A previous study investi¬ 
gated the effect of the wind termin ation shock assumin g 
a shock radius of 500-30,000 km (jKuroda et al.l[2?)0^ . 
They suggested that, in an ordinary SN, the wind ter¬ 
mination jumps appear only when the temperature suf¬ 
ficiently decreases, so the wind termination has little ef¬ 
fect on the final abundance pattern. However, the wind 
termination shock in our HN wind may appear at suf¬ 
ficiently high temperature because of the high neutrino 
luminosity from the mPNS. Thus, it can affect the final 
abundance pattern. 

In the neutron-rich wind, when the post-shock tem¬ 
perature is high enough (Tg > 2 x 10®K), the heavy 
elements can be photo-disintegrated into lighter nuclei. 
Therefore, the wind termination shock can hinder the r- 
process. In Fisrure fTHl we show the nucleosynthesis result 
of the trajectories of different shock radii (3000, 5000, 
7000, 10,000, and 15,000 km). The other parameters of 
the trajectories are the same as trajectory (a) (denoted 
as “w/o shock” in Figure fT6l). The abundances of the 
nuclei heavier than the second peak become smaller for 
smaller radius of the termination shock. This is because, 
for a shock radius smaller than 10,000 km, the post-shock 
temperature increases to be > 1.5 x 10®K and is sus¬ 
tained for a long time, and consequently, the elements 
larger than the second peak are photo-disintegrated into 
lighter nuclei. In the case where the shock radius is larger 
than about 10,000 km, on the other hand, the post-shock 
temperature is below 1.5x 10®K and the effect of the wind 
termination shock on the final abundance is small. 

5.4.2. The Effect on the vp-process 

The z/p-process continues longer in the presence of 
the wind termination because the temperature decreases 
more slowly compared to t he case withou t the w ind ter¬ 
mination. As explained in iWanaio et al.l (j2nilh . strong 
z/p-processes occur if the post-shock temperature be¬ 
comes Tg > 1.5 — 3. Such a high post-shock temperature 
may be achieved in our HN wind model because of the 
higher neutrino luminosity. Therefore, in such a case, the 
z/p-process may proceed more efficiently than the wind 
without the termination shock. From the viewpoint of 
heavy p-nuclei production, it is worth investigating the 
effect of the wind termination to study the synthesis of 
these elements in our HN wind model. 

In Figure 1171 we compare the final abundance for the 
trajectories with different shock radii. First of all, note 
that the abundance pattern for A < 80 is significantly 
affected by the presence of the shock termination because 
proton-capture reactions in the region of smaller atomic 
number can proceed for a long time owing to the higher 
temperature. 

For the heavier region, on the other hand, the effect 
of the shock termination depends on the location of the 
shock radius. When the shock radius is smaller than 
10,000 km, the post-shock temperature becomes high 
enough so that additional proton-capture reactions oc¬ 
cur as expected above, and consequently, the z/p-process 
can produce heavier elements. In addition, since the fluid 
elements are decelerated by the shock, they are exposed 


to large neutrino flux for a long time. As a result, the 
z^'p-process proceeds for a long time. Indeed, in the tra¬ 
jectories with the shock radius at small radii i?sh = 5000 
and 7000 km, the final abundance patterns are modified 
by the long-lasting r-p-process. In these trajectories, the 
heavy p-nuclei that are underproduced in the y-process 
are indeed synthesized. ^^®La is produced in the trajec¬ 
tory of i?sh = 5000 km and are synthesized in 

the case of i?sh = 7000 km. Thus, the HN wind model 
with the wind termination shock can be a possible pro¬ 
duction site of these p-nuclei. 

More interestingly, in the case of i?sh = 3000km, the 
final abundance pattern totally changes from the result 
of the trajectory without the termination shock. In this 
trajectory, the elements up to A ^ 150 are once pro¬ 
duced. However, the elements of A ^ 110 — 140 are 
photo-disintegrated because of the high post-shock tem¬ 
perature (Tg ~ 2). The elements of the neutron magic 
number N = 82 are relatively stable, so they avoid photo¬ 
disintegration. As a result, the final abundance pattern 
has a sharp peak around A = 145. 

Note that the final abundance pattern for the vp- 
process is less affected by the shock termination for the 
larger shock radius because the temperature range rel¬ 
evant to the z/p-process is higher that that for the r- 
process. 

5.5. The Effect of the Rotation and the Magnetic Field 
of the PNS 

In general, PNSs have an angular momentum and 
a magnetic field. If the PNS has a strong magnetic 
field and rotates rapidly, the magnetic force is impor¬ 
tant for the dynamics of the win d and for the ther¬ 
mal history of the flui d element (IM etzger e t al.l I2007t 
iVlasov et a.1.11201^ . In IMetzger et al.l (|2007ll . the crite¬ 
rion for which the wind is magnetically driven is shown 
(see their discussion above Equation (17) in their paper). 
Using our parameter set (i.e., 100, Ri, = 15 

km) and the quantities obtained our wind solutions (the 
asymptotic velocity v°' ~ 0.3c and the mass-loss rate 
M 4 X 10“^ Mq s“^), we can estimate the rotation 
period below which our HN wind is magnetically driven 
as 

P Ri 0.6 Bi 4 ms, (49) 

where B 14 is the magnetic field strength at the neutri- 
nosphere in units of 10^^ G. This is very close to the 
break-up period ^ 0.6 ms. Therefore, the effect of the 
rotation and the magnetic field on the wind is negligi¬ 
ble unless the rotation period is approximately the same 
as the break-up period or the magnetic field is strong 
enough. 

If the short rotational period and the strong magnetic 
field are achieved, the temperature-decrease timescale 
can become even smaller. (It can be below 1 ms depend¬ 
ing on the field strength and th e rotational period as de¬ 
scribed in IMetzger et al.1 (|2007[1 .) However, the fluid ele¬ 
ment accelerated by the magnetic force can pass through 
the neutrino-heating region quickly. Therefore, the mag¬ 
netically driven wind has less asymptotic entropy than 
that of the wind driven only by neutrinos. It is important 
to investigate which of these opposite effects is dominant 
in our HN wind, but this is beyond the scope of this 
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Figure 16. Abundance distribution of the elements as a func¬ 
tion of mass number A. The final abundance patterns of the 
trajectories of different shock radii (3000, 5000, 7000, 10,000, and 
15,000 km) are shown. The other parameters of the trajectories 
are the same as in trajectory (a), the result of which is also shown 
in the figure. 

work. We will investigate these effects next time. 

5.6. The Effect of Fallback 

The wind material could fall back to the central rem¬ 
nant (neutron star or black hole) with infalling stellar 
material. Therefore, the ejected mass of our HN wind 
can be somewhat smaller than that estimated in the pre¬ 
vious sections. 

A numerical study abo ut aspherical (je t-induced) ex¬ 
plosion of a Pop III star (iTominasral I2009D suggests that 
if the energy deposition rate is high enough, the mate¬ 
rial just above the central remnant can be ejected toward 
the jet axis. We consider the mass ejection toward the 
polar region so the wind material can avoid falling back 
because the ram pressure of the infalling stellar matter 
is expected to be weak in the region. 

6. CONCLUSIONS 

In this work, we investigated possible processes of 
heavy-element synthesis in HN explosions of very mas¬ 
sive (IOOMq) stars. We considered the neutrino-driven 
winds from the inPNSs formed after the collapse. We 
constructed a wind solution for the parameter set tabu¬ 
lated i n Table [J based on th e numerical relativity simu¬ 
lation (|Sekiguchi et al.l[2012D and calculated nucleosyn¬ 
thesis on the winds. We parameterized the average ener¬ 
gies of electron neutrinos and electron antineutrinos 
and investigated heavy-element synthesis processes. 
Here we make our conclusions as follows. 

1. The temperature-decrease timescale at T = 
0.5 MeV of the mPNS wind (about 6 — 7 ms) is 
much shorter than that of the ordinary PNS wind 
(about 40 — 50 ms), and thus the environment of 
mPNS wind is conducive to heavy element synthe¬ 
sis. 

2. In the neutron-rich region of the parameter space, 
although strong (i.e., making all elements up to 
third peak) r-process hardly occurs, weak r-process 
occurs and the second-peak elements are mainly 
synthesized. Near the fiducial parameter, there is 
a parameter region in which the wind reproduces 
the weak r elemental abundance pattern of some 
metal-poor stars, e.g., HD 122563. The third-peak 
elements can be produced only in the narrow pa¬ 
rameter region Cp^ — ^Pe ft MeV. 



mass number A 


Figure 17. Same figure as Figure 1161 but in the case of the 
proton-rich trajectories. The final abundance patterns of the tra¬ 
jectories of different shock radii (3000, 5000, 7000, 10,000, and 
15,000 km). Other parameters of the trajectories are the same as 
trajectory (e). 

3. In the moderately proton-rich region, i^p-process 
occurs and some light p-nuclei (e.g., proton-rich iso¬ 
topes of Ru, Pd, and Cd) are synthesized. In some 
parameter region around Cp^ ss Cp^, strong ( 71 , 7 ) 
reactions proceed and make the abundance distri¬ 
bution into a more neutron-rich one. As a result, 
the neutron-rich elements between A = 130 — 140 
are mainly produced. 

As discussed in the previous section, the uncertainties 
in the other parameters such as Lp., M, and Rp, which 
were unchanged in this work, are also expected to affect 
the nucleosynthesis in the wind from mPNSs. In the next 
study, we will discuss the possibility of HNe to the Galac¬ 
tic chemical evolution through the parameter survey of 
the nucleosynthesis in the wind about these parameters. 
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APPENDIX 

DETAILS OF THE REACTION AND COOLING RATES 


In this appendix we summarize the cooling and reaction rates used to construct the neutrino-driven wind solutions. 
The cooling rate due to electron and positron capture by free nucleons is given according to FFN85 as 


In 2 

QeN — F ‘^p^e\ec) i 


(Al) 


where rip and n„ are the number densities of protons and neutrons, respectively, and (/t) is the /t-value of the 
transition between nucleons, given as 

2-k^ 

(ft) = In 2 • „ „ -—^ 

Gl{gv + ig\)ml 

ft! 1.0157 X 10^ s (A2) 

where Gp is the Fermi coupling constant and gv = ^ an d qa = 1.275 are th e vector- and axial-vector-type coupling 
strength, respectively (the values of them are taken from [Piaster et alll2012ll . Je is given using phase-space factor J 
as 


Je=(—) r dr?, 77^(17, + C„)' 

\me J JrfL 


1 


1 - 


1 


fkBT 
V l^e 


1 


1 - exp(77^ - Cn - jyf) 


exp(r 7 , - r]f) + l [ exp(r 7 e -b Cn - + 1. 


J. 


(A3) 


The phase-space factor J is given as 

J= - Ve) - F 5 {r}u “ C« “ Ve) 
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+ l^Ve + Kn][F4{Ve “ ) “ Fi{r}u “ Cn “ 

+ [10(,7f')2 + 12,y,^Cn + 3(Cn)^][i^3(r;f - V^) - F^ - Cn - 
+ + ^Ve iCnf + (Cn)^] [^2 (^if “ Ve) “ F2iv^ “ C« “ ??e )] 

+ [5(,7i^)4 + 12(,7i^)3c„ + 9(,7i^)2(C„)^ + 2r?f'(C„)3][i^i(77f - - F^iv^ - (n - V^)] 

+ [iVef + KVe)‘^Cn + KVefiCnf + (Cn )^] [^0 (^if “ Ve) “ FoiVu “ Cn “ Ve% 

(A4) 

where 77 ^ is A/U-bT in the case of electron capture by protons and is rUe/kBT in the case of positron capture by 
neutrons, and Fk(rf) is the Fermi integral defined as 


POO 

Fkiv) = dx 
do 


1 


exp(a; — ry) -|- 1 


(AS) 


which is evaluated using the approximated expressions in FFN85. 

Th e cooling rate per unit m ass due to electron-positron pair annihilation into a neutrino-antineutrino pair is given 
from ICooperstein et al.l (|1986D as 

Qe'^e~ ~ ~F\oss^e+e~ : (A 6 ) 

P 


where E\oss is the average energy of the neutrino-antineutrino pair emitted per reaction and given as 

P l'^4(?7f) , A4(-77f)' 

-^loss ' TP / F\ ’ 

L^3(?7e) ^3(-?7e)J 

and Fg+g- is the annihilation rate of the electron-positron pair per unit volume and given as 

_ Gp(Dg -|- + Dr) - 


® ® 97r 

C/gT is the energy density of electrons and positrons given as 

(fcBT)4 


“ g— 


F3i±v[), 


(A7) 


(AS) 


(A9) 


and i3g+e- is the final state blocking factor given using the average energy of emitted neutrino E^, = £'ioss /2 as 

-1 r / f=, \ T -1 


Sg 


1 + exp 77 ^ - 


E, 


kBf) 


1 + exp 77 p - 


E, 


kBf) 


(AlO) 


This rate is derived using the approximation that the distributions of electrons and positrons are in thermal equilibrium 
(i.e., the chemical potential of positron is —y^g). Equation (IA 6 I) can be written as 


lG'^(Z?e-|-i:)^-|-£),-) g 

Qe+e- = - --(feT) K. 

where the phase-space factor K is given as 

K = [F4{Pe)F3{-Pe) + F4{-Pe)F3{r]e)\Be+e-- 

The reaction rate of electron and positron capture by free nucleons Aec, Ape is given from FFN85 as 

X-—I 

(/*)■’ 

where /g is given using the phase-space factor of this reaction I as 


(All) 

(A12) 

(A13) 
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1 - exp(? 7 ,f - Cn - vF 


(AM) 


and I becomes 


I = Fi{r]e - Pe) - Fiiv^ - Cn - Ve) 

+ (2C„ -k 477f')[F3(?7f - 77 ^) - F^irj^ - Cn - Ve)] 
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+ + ^CnVe + C]{P2{Ve “ Ve) “ F2{v^ “ Cn “ Ve)] 

+ + eiv^rCr. + 2v^iCnnF,{r,f - r,^) - - C„ - V^)] 

+ + ‘^CniVef + iCnf iVe fWoiVe “ Ve) “ Foilu - Cn - Ve)]- (A15) 


THE TREATMENT OE THE WIND TERMINATION SHOCK 


In this appendix we describe the treatment of the wind termination shock. Using the adiabatic sound speed = 
{dP/dp)s and the adiabatic gamma 7 = clKPjp), the specific enthalpy (including the rest mass) h and the pressure 
P can be expressed as 


h = l 


Cs p _ pet 
7-1 7 


We dehne the ratios x and y as 


Ps Ps 

x = — = —, y = -^, 

Pw Us Jw 


then Equations (ITtI) and 


become 

1 + 7wM.^( 1 - 1/a:) + 7wMw/(l ~ 7w) 

1 + (l/a;2)7sM2_ 1) 

+u'i)-x‘^ [1 


y = - 


0=1 + 2 / 


^S,W 


7w 7s - 1 


7w 1 


(1 + 


where M = ujcs is the Mach number of the fluid. These equations are combined to give the equation of x as 


7s 


7w 7s - 1 


Jv/M^ + 1 + 


7sMw 
7s - 1 


+ u^) - 1 + 


7w - 1 


(1 + I x^ + 


7sMw 
7s - 1 


= 0 . 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


This equation has four solutions for x. One of them gives the true post-shock quantities. Another solution corresponds 
to the unshocked (no discontinuity) solution because this solution becomes x = y = 1 if we assume 7 s = jw The 
other solutions are not real, so we do not consider them. The solution related to the wind without the shock is close 
to unity, so we can obtain the post-shock solution by finding the real solution of Ea. (IB5|) greater than another real 
solution, i.e., the unshocked one. 

In this method, we assume the adiabatic gamma of the post-shock material 7 s. In general, it does not coincide 
with the adiabatic gamma evaluated fr om t he obtained post-shock quantities. In order to find the true post-shock 
quantities, we iteratively solve Equation (|B5p until the assumed % coincides with the one evaluated from the post-shock 
quantities. 












